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In nonlocal general relativity, linearized gravitational waves are damped as they 
propagate from the source to the receiver in the Minkowski vacuum. Nonlocal grav- 
ity is a generalization of Einstein's theory of gravitation in which nonlocality is due 
to the gravitational memory of past events. That nonlocal gravity is dissipative is 
demonstrated in this paper within certain approximation schemes. The gravitational 
memory drag leads to the decay of the amplitude of gravitational waves given by the 



O 

exponential damping factor exp (— t/r), where r depends on the kernel of nonlocal 

U 

gravity. The damping time r is estimated for gravitational waves of current obser- 
vational interest and is found to be of the order of, or longer than, the age of the 
Q\ ■ universe. 
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I. INTRODUCTION 

A classical nonlocal generalization of Einstein's theory of gravitation |l| has been devel- 
oped in recent papers [2h3] . In this theory, nonlocality is introduced into gravitation theory 
via a causal scalar "constitutive" kernel that acts as the weight function for a certain av- 
erage of the gravitational field over past events. The gravitational field is represented by a 
local field that satisfies integro-differential field equations. Thus classical nonlocal gravity 
is Einsteinian general relativity endowed with a certain simple memory of past events; this 
nonlocal aspect of gravity simulates dark matter. That is, in this theory, nonlocality is the 
main source of phenomena associated with what appears as dark matter in astrophysics. 
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To account for the "flat" rotation curves of spiral galaxies [8HlO|. the characteristic length 
scale of nonlocality is assumed to be of order A = 10 kpc. This is also the scale over which 
memory of distant events fades. 

It follows from the weak field approximation of general relativity that gravitational waves 
satisfy the ordinary wave equation and thus propagate in vacuum with the speed of light. 
However, in the treatment of linearized gravitational waves in nonlocal general relativity, 
one encounters the following nonlocal wave equation for the free propagation of linearized 
gravitational radiation [7J 

Dhij + Sij(x) + / R(x - y)s ij (y)d i y = , (1) 

where hy is the amplitude of gravitational radiation and 

s ij {x) = -J K(|x-y|,x-y)^(ci-|x-y|,y)tfy (2) 

Here x represents an event in spacetime with coordinates x M = (ct,x) and □ := "i] a ^d a dp is 
the d'Alembertian (wave) operator. Greek indices run from to 3, while Latin indices run 
from 1 to 3. Moreover, g^ u = r]^ u + h^ u is the metric tensor of spacetime and the Minkowski 
metric tensor r] a p is given by diag(— 1, 1, 1, 1) in our convention. We use units such that 
c = 1, unless otherwise specified. The freedom in the choice of coordinates, which translates 
into the gauge freedom of dimensionless gravitational potentials h^, can be employed to set 
hofj, = 0, impose the transverse gauge condition dh^/dx^ = and render hij traceless |7|. 
The gravitational wave amplitude h^ thus characterizes the deviation of the spatial metric 
away from flat 3D Euclidean space. 

It should be emphasized that our way of introducing nonlocality into general relativity is 
by no means unique; however, it is the simplest approach to nonlocal general relativity that 
is widely consistent with observation [2nZ|- 

In the general linear approximation of nonlocal gravity, the scalar constitutive kernel of 
the theory is a universal function given by the causal convolution kernel K(x — y). Its 
reciprocal R(x — y) is a causal convolution kernel as well, as discussed in detail in Ref . [7| . 
The reciprocity between K and R implies that 

K(x -y) + R(x -y)+ K(x - z)R{z - y)d 4 z = . (3) 

It is interesting to note that K and R can be interchanged in the integrand of Eq. ([3]) by 
simply changing the variable z to z' given by z' = x + y — z. We find it useful to work in the 
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Fourier domain, which is permissible for the functions under consideration here [6|. Thus 
let 

/(0 = J ' f{x)e-* x d A x (4) 

be the Fourier transform of /, where £ • x :— r] a /3C, a x' 3 . Then, 

/(^) = ^l//(Oe^'A. (5) 

It follows from Eq. (J3J) and the convolution theorem that in the 4D Fourier domain, 

K + R+KR = 0. (6) 

Given K, one can therefore determine R using Fourier transforms and vice versa. These 
kernels have dimensions of (length) -4 , so that K and R are dimensionless. 

Q 

For the reciprocal kernel R, we adopt the expression derived in Ref. [7j], namely, 

R(x -y)= H(x° - y° - |x - y\)A e - A ^%(x - y) . (7) 

Here H(s) is the Heaviside unit step function such that H(s) = 1 for s > and H(s) = 
for s < 0, 1/A is a constant length such that < AXq < 1, where Ao = 10 kpc is the 
basic nonlocality scale in this theory, and q > is the reciprocal Newtonian kernel. This 
reciprocal kernel occurs in the Poisson equation for the gravitational potential V of nonlocal 
gravity in the Newtonian regime, namely, 



V V(x) = 4ttG [p(x) + J g(x - y)p(y)rf 3 y] . (8) 



Thus V is due to matter of density p and "dark matter" of density pr>, which is the con- 
volution of p with the reciprocal Newtonian kernel q; that is, nonlocality simulates dark 
matter J8l-ll0|. The reciprocal kernel q has been discussed extensively in Ref. |6(, where two 
possible examples were worked out in detail, namely, 

v^ e~ Ar , (9) 
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^ = ir T,r e '" (10) 

47tA r(a + r) 
Here r = |r| and a, a/Ao <C 1, is a constant length. The kernels q' and q" are real positive 

functions that are integrable as well as square integrable over all space; moreover, we will 

assume here that their Fourier transforms are real positive functions as well. This was proved 
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for q" in Ref. [6|, but the situation regarding q' is more complicated. Further discussion of 
this issue is contained in section V, where it is argued that q' is positive for sufficiently small 
values of a/Ao- In fact, it follows from the numerical results presented in figure 1 of Ref. 6J 
that q' > for a/Ao = 0.001. Therefore, we will assume in this paper that a/Ao is always so 
small compared to unity that q' is positive. 

The values assigned in Ref. |6| to the constant parameters Ao, A and a of kernel (J7|), 
namely, Ao = 10 kpc, A\ = 0.1 and a/Ao — 0.001, are rather tentative and will be used 
in this paper for some of the numerical estimates in section V. In general, Ao should be a 
galactic length scale to account for the rotation curves of spiral galaxies, while A\ Q and a/Ao 
should be positive and small compared to unity such that < a A <C 1. 

Let us note here that R is a positive function that is integrable as well as square integrable 
and satisfies the causality requirement, namely, it is nonzero only when x^ — y^ is a future 
directed timelike or null vector in Minkowski spacetime. This means that event y must be 
within or on the past light cone of event x such that x° > y° + |x — y|, where |x — y| can 
be characterized as retardation. The calculation of R becomes much simplified if we neglect 
retardation; then, 

%k)«/-g(k), (11) 

where we henceforth use the approximation sign instead of the equality sign to indicate that 



the calculation neglects retardation and is therefore not exact. Ref. 



7j should be consulted 



for a detailed discussion of the approximation scheme that is basec 
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on neglecting retardation. 



Using Eq. (Q and inverse Fourier transformation, we find that [7] 

K (x°, x) « --^—H(x°) I q(k)e ik *e- A ( 1+ V x ° d 3 k . (12) 

(2-rcy J 

The choice of a particular gauge for the gravitational potentials together with the require- 
ment of causality implies that the explicit form of our nonlocal gravitational wave equation 
is valid in a particular inertial frame and violates time-reversal invariance. 

To understand the physical import of Eqs. (OQ) and (J2J), a mechanical analogy turns out 
to be quite useful: These equations are reminiscent of the equation of motion of a linear 
oscillator with a dissipation term that is proportional to the velocity of the oscillator. Here 
dhij/dt is suggestive of the "velocity" of the oscillator. With the appropriate sign for the 
coefficient of the dissipation term, one has a damped oscillator. In a similar way, with the 
proper functional forms for the nonlocal kernels, Eqs. ([T]) and (J2J) indicate free propagation 



of gravitational waves with damping. The nonlocality of the theory originates from a certain 
average over past events in spacetime; this memory of the past thus appears to act as a drag 
that dampens the free propagation of linearized gravitational waves. In nonlocal gravity, 
memory fades exponentially for events that are distant in space and time; see Eqs. (J7|), ([9]) 
and (II Oft . Similarly, the amplitude of gravitational radiation decays exponentially with time 
as exp (— t/r), where the damping time r is related to the nonlocal kernel. 

In general, as the waves propagate freely through Minkowski vacuum, the wave amplitude 
may grow or decay in time due to nonlocality. We expect that, with the correct nonlocal 
kernel, the solutions of our linear nonlocal homogeneous wave equation are well behaved and 
decay in time leading to the stability of Minkowski spacetime under small perturbations. 
We show in sections II and III that there is no instability and the waves indeed decay in 
time when retardation is neglected in the kernel; that is, employing the same approximation 
scheme as in Ref. (7| and using the proper physically reasonable nonlocal kernel, we find 
that the waves are exponentially damped. A simplified 2D toy model is discussed in section 
IV, where exponential damping is demonstrated even in the presence of retardation. It thus 
appears from the physical arguments provided in this paper that damping of gravitational 
waves is a feature of nonlocal gravity; however, the general mathematical problem involving 
retardation remains unsolved. The damping time r is estimated in section V for the case 
of gravitational waves that are of current observational interest. Section VI contains a brief 
discussion of our results. 

II. NONLOCAL WAVE EQUATION 

Let us now concentrate on a component of the gravitational wave amplitude and replace 
Eqs. (II]) and (J2J) by the equivalent set 

Dip + S(x) + I R{x- y)S(y)d A y = (13) 

S(x) = — J A'(|x - y |, x - y) if)(t - |x - y |, y) d 3 y . (14) 

As these equations are linear in ip and the kernels are real, it is convenient to work with 
a complex wave amplitude if) with the understanding that the real part of if) is physically 



and 



significant. We then look for a solution of the form 

*{?) = e" Ja V(x) , (15) 

where a is in general complex. Substituting our ansatz in Eqs. ( 1T31) and ( Tl4l) . we find 

(V 2 + « 2 )0(x) +ia J W a (x - y)0(y) d 3 y = , (16) 

where W a is given by 

W a (z) := F a (z) + — — / e (--A)| z - u | g(z _ u ) Fa(u)rf 3 u (17) 

/i — la J 

and 

F Q (z):=-e 4Q l z lir(|z|,z). (18) 

The main issue here is whether all solutions of Eqs. (Ti~6l)-(TT8l) that satisfy proper boundary 
conditions are such that a, 

a = u + iA, (19) 

has a negative imaginary part A < 0. In this case, 

^( x ) = e -*"V(x) e At (20) 

will exponentially decay in time. Otherwise, the perturbation will blow up as t — > oo, 
which is physically unacceptable, as it would indicate an intrinsic instability of Minkowski 
spacetime within the framework of nonlocal gravity. 

In the Fourier domain, if 0(k) is nonzero, Eq. ( 1T6|) can be written as 

a 2 - \k\ 2 + iaW a (k) = 0, (21) 

which expresses the dispersion of gravitational waves due to nonlocality. Writing a as in 
Eq. f|T9|) and introducing the real (271) and imaginary (2X) parts of W a , 

W a = 2(71 + it) , (22) 

we find that Eq. ( 12 ip reduces to its real and imaginary components 

w 2 -A 2 -|k| 2 -2(u;Z + Aft) = (23) 

and 

wA + (cuTZ - AX) = . (24) 



It follows from Eq. ((24} that 

A -~^x' < 25 > 

which can be substituted in Eq. ( )23|) . In the resulting equation, we can choose u — X as a 
new variable and after some algebra we find 

u = 1 ± ( J + VJ 2 + ^ 2 ^ 2 ) , (26) 

where J7" is given by 

.7 = -(|k| 2 -^ 2 +X 2 ). (27) 

If A turns out to be negative independently of the sign of u, then the solution of the nonlocal 
wave equation will decay in time and could be physically acceptable. As discussed in the 
following section, thus far it has been possible to show this for Eqs. (fT6|) - (TT8|) only in the 
approximation that retardation is neglected. 

III. DAMPING 

Consider the approximation scheme, introduced in Ref. [7[, that involves neglecting the 
retardation in Eqs. ( IT71) and ( {TBI) . This means that Eq. ( IT8l) is replaced by F a (z) ps —K(0, z). 
Taking Eq. ( 1T21) into account, we find that in this approximation F a (z) is independent of a, 
since 

F a (z)^Aq(z). (28) 

Similarly, Eq. f|T7|) is approximated by 

W a (z) « Aq{z) + -A— / q(z - u)q(u)d 3 u , (29) 

A — tot J 

so that its Fourier transform can be expressed as 

W a (k)^Aq(k)\l + 4^-}- (30) 

L A — ia\ 

We recall that A is a positive parameter, < AA < 1, and g(k) > by assumption. It 
is shown in section V that g(k) is a function of |k| and is real and positive when a/Ao is 
sufficiently small compared to unity. It follows from Eq. ( 1221) that 



n~ l -Aq 



A(A + A) ; 
(A + A) 2 + u 2 q 



(31) 



and 

Substituting these results in Eq. ( |25|) . we find after some algebra that A can be expressed 

as 



1 a* 
Ak — Aq 

2 H 



A 1 

1+ , , , A xo , — o q 



(A + Ay + u^r (33) 

which means that A < and all the modes decay regardless of the value of u. Indeed, 
Eq. (133]) can be expressed as a cubic equation for A with all positive coefficients, since 
A > and q > 0; therefore, there is at least one real root, which must be negative, as it 
follows from Descartes' rule of signs that the cubic equation cannot have any positive root 
for A. 

We remark here for the sake of completeness that the substitution of Eqs. ( |3T1) and ( 132]) 
in Eq. ( 126]) results in a second formula involving uj and A. This equation and Eq. (133]) then 
constitute two coupled algebraic equations for the two unknowns u and A; in principle, one 
can determine u and A in this way in terms of |k|, g(k) and A. 

To what extent does the result that A is negative depend upon neglecting retardation? 
It is shown within the framework of a simple 2D toy model in the next section that all the 
modes decay even in the presence of retardation. 

IV. TOY MODEL 

Consider the 2D nonlocal wave equation given by 

n^-e— x(x-y)V{t-\x-y\,y)dy = 0, (34) 

where e, < e < 1, is a small parameter. Working to first order in e, we wish to show 
that all the modes decay provided kernel x( z ) is the Fourier transform of an even positive 
function x(C), 

X(C) = X(-C)>0, (35) 

which implies, among other things, that the function x( z ) * s rea l an d even. Therefore, 

X(C) = 2 / X(z) cos ((z) dz , X {z) = ~ / X(0 cos (C*) d( . (36) 

./n TT In 
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The assumptions regarding kernel x are reminiscent of the fact that in nonlocal gravity, 
kernel q, for instance, is such that q > is just a function of |k|; in fact, kernel x is the 2D 
analog of Aq in Eq. ( 1281) of the previous section. 
We assume, as before, that 

V{t,x) =e" iaf $(x), (37) 

so that Eq. (1341) takes the form 

d 2 &(x) f 

\ ' + a 2 $(x) + iea / I(x - y)$(y) dy = , (38) 

where 

I{z) := e^x(z) ■ (39) 

Nonlocal wave equations of the general type of Eq. ( |55|) have been the subject of previous 



investigations; see, for instance, Ref. 11] and the references cited therein. 

Let $(/c) be the Fourier transform of $(x); then, it follows from Eq. (1381) that if $(fc) is 
nonzero, we have 

a 2 - k 2 + ieal{a, k) = , (40) 

where 

/(«,*) = /**•*<-.«<,. (41) 

Working to first order in e in the Fourier domain, it is clear from the treatment of section 
II that for a = u + % A, we have 

w = ±fc+-e/j(a ,fc), A = --e/ fl (a ,fc), (42) 

where Jj and J/j are the imaginary and real parts of /, respectively, and «o := ±fc is the 
solution of Eq. fj4*0l for e = 0. To prove that to first order in e all the modes decay, we need 
to show that Ir{ocq, k) > 0. 

To this end, we first note that the real part of Eq. (T4T1) for a = a involves cos (a \z\ — kz), 
since x( z ) is rea l; hence, 

x(z) cos («o + k)z dz + / x( z ) cos ( a o — k)z dz . (43) 

oo Jo 

Changing the variable z to — z in the first integral, we find 

POO 

In(a ,k)= / x( z ) cos (a + k)z + cos (a — k)z dz . (44) 

Jo L -I 
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Substituting either k or — k for a in Eq. (144]) . the integrand remains the same; hence, 

POO 

In{otQ,k) = \ (1 + cos2kz)x(z) dz . 
Jo 



It follows from this result and Eq. (136]) that 



I R (a ,k) = - x(0) + xC2k) 



But x is positive by assumption, hence 



A = --e[x(0) + x(2fc) 



<0 



(45) 



(46) 



(47) 



which means that to first order in e all the modes indeed decay even in the presence of 
retardation. 

We now return to Eq. ( |33l) and use this result to estimate the nonlocality-induced damping 
time for gravitational radiation. 



V. DAMPING TIME r 



It follows from our approximate treatment in Eq. (|33|) that the damping time r := — 1/A 
depends upon q(k) as well as the frequency of radiation uj. It is intuitively clear that in 
nonlocal gravity, the propagation of gravitational waves with wavelengths comparable to, or 
longer than, the basic length scale Ao = 10 kpc could be significantly affected by nonlocality. 
In particular, a rough estimate for r involving gravitational radiation of wavelength compa- 
rable to the nonlocality length scale Ao would be a damping time of order Ao/c «3x 10 4 
yr. On the other hand, various observational efforts are under way to detect gravitational 
waves with wavelengths that are much shorter than Ao 12]. It turns out that in this case 
uj « ±|k| and that nonlocality generates only a small perturbation on the propagation of 
such waves J2|. In particular, current observational possibilities involve gravitational waves 
in the frequency range v > 10~ 8 Hz, where 2ttu ps |k|; therefore, Ao|k| > 6 x 10 4 for 
wave vectors of current experimental interest [l2J. For such radiation, A/|k| < 10~ 5 , since 
A\q < 1, and to calculate r from Eq. (|33l) we need a proper estimate for q(k). In fact, for 
the cases of observational interest, q(k) is small compared to unity. Moreover, the quantity 
in square brackets in Eq. (1331) is less that 1 + (A/u) 2 q, which is thus very nearly equal to 
unity. Therefore, it follows from Eq. ( 1331) that 

2 



Aq(k) 



(48) 
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for the frequency range that is the focus of observational searches for gravitational waves at 
the present time. 



The Newtonian regime of nonloca 



gravity is similar to the phenomenological Tohline- 



Kuhn approach to modified gravity 13M15I]; in particular, the reciprocal kernel g(r) is a 
generalization of the Kuhn kernel Q(t), 

w-ikh- (49) 

In fact, two examples of g(r), namely, q' and q" given respectively by Eqs. (jUJ) and ffTUl) . have 
been constructed starting from Q(t) in Ref. 6| by introducing positive constants a and A 
to change the functional form of the Kuhn kernel for r — > and r — > oo, respectively. It has 



[Q 



been shown in Ref. [6] via general arguments that the Fourier transforms of q' and q" are 
actually dimensionless functions of |k| and have the following properties: q" > 0, q" > q' 
and q' > — a/Ao; moreover, q' and q" both vanish as |k| —$■ oo. On the other hand, from the 
expressions for q' and q" given respectively by Eqs. (Q and (TTU1) . it is possible to compute 
explicitly the Fourier transforms of these functions using 

47T f°° 

g(k) — — / rq(r) sin(|k|r)dr . (50) 

|k| Jo 

Let us first consider the special limiting case of a = 0; then, q' and q" coincide and the 

resulting reciprocal Newtonian kernel qo, 

l 1 1 . M _ Ar 



has the same short distance behavior as the Kuhn kernel. It is straightforward to compute 
the Fourier transforms of Q and go using Eq. (|5T7|) and the results are 

<i^\ (52) 

and 

* = A (^ + W) + W\ aiCtan ^ • (53) 

Here we have used the fact that for real values of the constants p\ and p2, 

/ e~ pix sin (p 2X ) — = arctan (^) , Pl > ; (54) 

Jo x Pi 



see formulas 3.941 on page 489 of Ref. [161 ] . Thus go > 0, qo(0) = 2/(AXq), while for 



Aojkl ^> 1, g ~ Q- The similarity between the behaviors of go and Q for large wave 
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numbers in the Fourier domain is naturally related to the fact that go an d Q have much the 
same behavior as r — > 0. 

It should be emphasized that for a = 0, the Fourier transform of the reciprocal kernel is a 
positive function given by Eq. ( 1531) and this positive character of q is expected to persist for 
sufficiently small a/Ao <C 1. Indeed, we suppose that a/Ao is always so small that <f (k) > 0; 
rown numerically in figure 1 of Ref. [6| for a/Ao = 10~ 3 . Moreover, figures 
6( demonstrate that for AX = 10 _1 and a/\ = 10~ 3 , <f and q" have 
rather similar functional forms: they both start from finite positive values at |k| =0 and 
monotonically decrease to zero as |k| — > oo. 

For a > 0, we find from Eqs. ([9]), ( TTOl and ( |50l) via an integration by parts that 



in fact, this was s 
1 and 3 of Ref, 



and 



Here 



q =q 



A 



Ao 



;r 



e z E 1 {Z) 



1 



A (A2 + |k|2) A |k 



e z E 1 (Z) 



Z = a(A + i\k\) 



(55) 



(56) 



(57) 



9ft and 3 indicate real and imaginary parts, respectively, and E\ (z) is the exponential integral 
function defined for a complex number z with positive real part, 3ftz > 0, as |l7| 

„-zt 



E 1 (z) 



t 



dt. 



(58) 



Let us note that Z has a small positive real part a A, < a A <^C 1, and a positive imaginary 



part a|k| for k ^ 0. If k = 0, Z is real and 3 e z E\(Z) 



using the lemma connected with Eq. (29) of Ref. p] — that for complex Z, Q' e E\(Z) 
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0; otherwise, one can show- 



<0. 



Moreover, it follows from q" > q' and Eq. fl55|) that 3ft e z Ei(Z) > 0. Therefore, <f > 
whenever 



q" 



0<-<- 

A o 3ft e z E l {Z) 

We assume that a/Ao is always so small compared to unity that q' > 0. 
It follows from Eq. ( 158]) that in general 



(59) 



\e z E x {Z)\ <e aA E x (aA) 



(60) 
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where equality holds for k = 0. Moreover, using the expansion 17] 



E 1 (z) = -C-lnz-J2 { - ZR ^, (61) 

n=l 

where C = 0.577... is Euler's constant, we find that for < aA <C 1, 

E l (aA)^-C-\n(aA). (62) 

We therefore conclude from Eqs. fl55l) and fl56l) that in general 

0<q"-q' <^— Si(oA) (63) 

Ao 

and 

0< ^ Mgm + ^r' (64) 

It turns out that for wavelengths of current observational interest, namely, Ao|k| > 6 x 10 4 , 
the imaginary part of Z is much larger than the real part of Z\ in fact, QZ/^R.Z = yl _1 |k| > 
6 x 10 4 , since A~ x > X by assumption. It therefore follows from Eqs. f)63p and f)64p that for 
gravitational radiation with A |k| > 6 x 10 4 , the damping time ( jjgj) is always longer than 

2A |k| 
T ° " Ae^ El {aA) ' (65) 

For AAo = 0.1 and a/Ao = 0.001, we can use Eq. (jB"2"|) with In 10 ~ 2.3 to conclude that 
ro > 5 x 10 9 yr, which is about one third of the current estimate for the age of the universe. 
It is interesting that To is simply proportional to the frequency of gravitational waves. In 
fact, nonlocality-induced damping could become significant for cosmo logical gravitational 
waves with very low frequencies. Inspection of Eq. ( 165]) reveals that for realistic values of 
the parameters a and A, r is longer than the age of the universe for gravitational waves 
that might be detectable in the foreseeable future with laser interferometers. We recall 
that for such devices, the waves should have dominant frequency > 1 Hz for Earth-based 



and > 10 4 Hz for space-based antennas. However, detection of gravitationa 



waves with 



dominant frequency of several nHz may be possible with pulsar timing arrays [12|. Therefore, 
in connection with waves of low frequency v ~ 10~ 8 Hz, we will next examine the limiting 
cases where < a|k| <^C 1 and a|k| 3> 1, respectively. 

Consider first the case where ^sZ = o|k| ^C 1; an example is given by v ~ 10 -8 Hz 
and a/Ao = 10~ 8 . It follows that \Z\ ^C 1 and to compute Ei(Z) in this case one can use 
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expansion (I6TJ) . A detailed examination reveals that in this case 



which is much smaller than unity. Eq. (166|) is reminiscent of the expression (1521) for Q; that 
is, a/Ao is in this case so small compared to unity as to be essentially negligible. Hence it 
follows from Eq. ( HH1) that in this case 

r*—. (67) 

Thus for Ao|k| >6x 10 4 and AXq = 0.1, we find that r > 2 x 10 10 yr, which is nearly 1.5 
times the age of the universe. 

Next, we consider the case where QZ = a|k| ^> 1. An example is provided by the choice of 
parameter a in Ref. [6], namely, a/X = 10~ 3 , so that QZ > 60. Hence, we need to compute 
q for large \Z\ 3> 1. It is possible to develop an asymptotic expansion for Eq. ( 1581) by 
using e~ zt dt = —z~ l d(e~ zt ) and repeated integrations by parts. The resulting (divergent) 



asymptotic expansion for \z\ ^> 1 is given by 17 1 



n=0 

It thus follows from Eqs. ( )55l) and ( |56l) that for |Z| ^> 1, the dominant terms in q' and g" 
are much smaller than unity, and are of the form 

^a^W' r ~v^- (69) 

Using these asymptotic estimates, we find from Eq. (|48l) that the corresponding damping 

times would be 

, A a 3 |k| 4 /; 2A a|k| 2 

T ~^^-' T ^^4— • (70) 

For AX = 0.1, a/A = 0.001 and A |k| >6x 10 4 , one can estimate that r' > 2 x 10 15 yr and 

t" > 2 x 10 12 yr. Thus in each of the cases considered here, the nonlocality-induced damping 

of gravitational waves of current observational interest is insignificant as the corresponding 

damping time is of the order of, or longer than, the age of the universe. 

VI. DISCUSSION 

This paper is devoted to the study of an interesting feature of nonlocal gravity that shows 
up in the nonlocal wave equation for the propagation of free linearized gravitational waves in 
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a background global inertial frame of reference. The waves are damped as they propagate in 
vacuum. The corresponding nonlocal wave equation is reminiscent of a harmonic oscillator 
with a linear nonlocal damping term. In the case of the simple damped harmonic oscillator, 
damping or antidamping would depend on the sign of the damping coefficient. A similar 
situation holds in the nonlocal case, where the nonlocal kernel of the theory must be such 
that the gravitational waves are all damped and Minkowski spacetime is thus stable. Using 
certain simplifying assumptions and restricting parameter a, < a/Ao <C 1, of the kernel to 
be small enough such that the spatial Fourier transform of the reciprocal Newtonian kernel 
q is positive, we have shown that all the modes do indeed decay. However, for the linearized 
gravitational waves that may be detectable in the foreseeable future, the amount of damping 
would be negligible, as the damping time would be of the order of, or longer than, the age 
of the universe. 

Specifically, the nonlocality-induced damping time r is estimated in this paper for radi- 
ation of frequency v > 10 -8 Hz, which is the frequency range that is the focus of current 



observational efforts 12(. We recall that v ^> uq, where, vq = c/Aq « 10 12 Hz and Aq = 10 



kpc is the basic galactic length scale associated with nonlocal gravity. It has been shown 
in a recent work 7| that for v > 10~ 8 Hz, the nonlocal deviations from standard general 
relativity are negligibly small, a circumstance that is consistent with our estimate that the 
corresponding r is of the order of, or longer than, the age of the universe. On the other hand, 
linearized gravitational waves with very low frequencies, v < b>o, would be highly damped 
in nonlocal general relativity. 

Finally, it should be emphasized that all of our estimates for the damping time r de- 
pend, of course, on our assumptions regarding the functional form of the kernel of nonlocal 
gravity. Ultimately, the kernel of the theory must be determined through comparison with 
observational data. Therefore, our estimates for r may have to be revised as a consequence 
of the confrontation of nonlocal general relativity with observation. 
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